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Abstract. The error on a real quantity Y due to the graduation of the measuring 
instrument may be asymptotically represented, when the graduation is regular and 
fines down, by a Dirichlet form on R whose square field operator does not depend 
on the probability law of Y as soon as this law possesses a continuous density. This 
feature is related to the "arbitrary functions principle" (Poincare, Hopf). We give 
extensions of this property to M. d and to the Wiener space for some approximations 
of the Brownian motion. This gives new approximations of the Ornstein-Uhlenbeck 
gradient. These results apply to the discretization of some stochastic differential 
equations encountered in mechanics. 

Key words : arbitrary functions, Dirichlet forms, Euler scheme, Girsanov theorem, 
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stochastic differential equation. 

Introduction. 

The approximation of a random variable Y by an other one Y n yields most often a 
Dirichlet form. The framework is general, cf. Bouleau [3] whose results are recalled 
§1.1 below. 

Usually, when this Dirichlet form exists and does not vanish, the conditional law 
of Y n given Y — y is not reduced to a Dirac mass, and the variance of this conditional 
law yields the square field operator T. On the other hand when the approximation 
is deterministic, i.e. when Y n is a function of Y say Y n = r] n (Y), then most often the 
symmetric bias operator A and the Dirichlet form vanish, cf. Bouleau [3] examples 
2.1 to 2.9 and remark 5. 

Nevertheless, there are cases where the conditional law of Y n given Y is a Dirac 
mass, i.e. Y n is a deterministic function of Y, and where the approximation of Y by 
Y n yields even so a non zero Dirichlet form on L 2 (Py). 

This phenomenon is interesting, insofar as randomness (here the Dirichlet form) 
is generated by a deterministic device. In its simplest form, the phenomenon appears 
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precisely when a quantity is measured by a graduated instrument to the nearest 
graduation when looking at the asymptotic limits as the graduation fines down. 

The first part of this article is devoted to functional analytical tools that we need 
afterwards. We first recall the properties of the bias operators and the Dirichlet 
form associated with an approximation. Next we prove a Girsanov-type theorem 
for Dirichlet forms which has its own interest, i.e. an answer to the question of an 
absolutely continuous change of measure for Dirichlet forms. At last we recall some 
simple properties of Rajchman measures. 

The second part is devoted to the case of a real or finite dimensional quantity 
measured with equidistant graduations. The mathematical argument here is basi- 
cally the arbitrary functions method about which we give a short historical comment. 

Several infinite dimensional extensions of the arbitrary functions principle are 
studied in the third part. The first one is about approximations of continuous mar- 
tingales whose brackets are Rajchman measures. Then we consider the case of the 
Wiener space on which the preceding results may be improved and other asymptotic 
properties are obtained concerning the approximation of the Ornstein-Uhlenbeck 
gradient. Eventually we apply these results to the approximation of stochastic dif- 
ferential equations encountered in mechanics and solved by the Euler scheme. 

I. Functional analytical tools. 

1.1. Approximation, Dirichlet forms and bias operators. 

Our study uses the theoretical framework concerning the bias operators and 
the Dirichlet form generated by an approximation proposed in Bouleau [3]. We 
recall here the definitions and main results for the convenience of the reader. Here, 
considered Dirichlet forms are always symmetric. 

Let Y be a random variable defined on (fl, A, P) with values in a measurable 
space {E,~F) and let Y n be approximations also defined on (ft, A, P) with values 
in (EjJ 7 ). We consider an algebra T> of bounded functions from E into 1 or C 
containing the constants and dense in L 2 (E, J 7 , Py) and a sequence a n of positive 
numbers. With V and (a n ) we consider the four following assumptions defining the 
four bias operators 

f Vy? G £>, there exists A[<p] G L 2 (E, J 7 , P y ) s.t. V% G V 
\ lim n _ 0O a n E[( ¥ »(y n ) - <p(Y)) X (Y)] = Ey ["%>]*]. 

f G T), there exists A[ip] G L 2 (E, F,W Y ) s.t. Vx G V 
\ lim n ^ 00 a n E[( V 9(r) - <p(Y n )) X (Y n )} = E Y [A[<p] X }. 

Vy? G V, there exists A[tp\ G L 2 (E, J 7 , P y ) s.t. Vx eV 
lim n ^ 0O a n E[(^(y n ) - <p{Y)){ X {Y n ) - X (Y))] = -2E Y [A[<p] X ]. 

Vy? G V, there exists A[(p] G L 2 (E, J 7 , P y ) s.t. Vx G V 
lim^oo aJZ[(<p(Y n ) - <p(Y))( X (Y n ) + X (Y))] = 2E Y ft[<p] X ]. 
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We first note that as soon as two of hypotheses (HI) (H2) (H3) (H4) are fulfilled 
(with the same algebra V and the same sequence a n ), the other two follow thanks 
to the relations _ _ 

7 A + A vi A ~ A 

A = = A = . 

2 * 2 

When defined, the operator A which considers the asymptotic error from the point 
of view of the limit model, will be called the theoretical bias operator. 

The operator A which considers the asymptotic error from the point of view of 
the approximating model will be called the practical bias operator. 

Because of the property 

( a [v},x)l2(f y ) = (ip,A[x})l*(v y ) 

the operator A will be called the symmetric bias operator. 

The operator A which is often (see theorem 2 below) a first order operator will 
be called the singular bias operator. 

Theorem 1. Under the hypothesis (H3) 7 

a) the limit 

~ ct 

8[tp, X ] = lim -?E[(tp(Y n ) - ip(Y))(x(Yn) - x(Y)] (1) 

defines a closable positive bilinear form whose smallest closed extension is denoted 
(£,©). 

b) (£,©) is a Dirichlet form 

c) (S, D) admits a square field operator T satisfying Vy?, x^V 

I>] = A\^\ - 2yA[y] (2) 
E Y [T[<p]x] = ]ima n E[(<p(Y n ) - V {Y)f{ X {Y n ) + X (Y))/2] (3) 

n 

d) (£, D) is local if and only ifW^p G V 

\ima n E[^(Y n )-^(Y)) 4 ] =0 (4) 

n 

this condition is equivalent to 3A > 2 lim ra a n E[|</?(l^) — </?(y)| A ] = 0. 

e) If the form (£, B) is local, then the principle of asymptotic error calculus is 
valid on V = {F(f f p ) : fi E T>, F e C 1 (1R P , R)} i.e. 

hm„ a n E[(F(/!(y n ), • • • , f P (Y n )) - F(f 1 (Y), f p {Y)f] 

= Ey e&=i FKh, ■ ■ ■ , f P )F;(h, . . . , u)m, m. 

An operator B from V into L 2 (Py) will be said to be a first order operator if it 
satisfies 
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Theorem 2. Under (HI) to (H4). If there is a real number p > 1 s.t 
lima n E[(^(y n ) - <p{Y)) 2 \1>(Y n ) - ^(Y)\"] =0 V</?, ip E V 

n 

then A «s /irs£ order. 

In particular, if the Dirichlet form is local, by the d) of theorem 1, the operator 
A is first order. 

1.2. Girsanov-type theorem for Dirichlet forms. 

An error structure is a probability space A, P) equipped with a local Dirichlet 
form with domain D dense in L 2 (Q,A, P) admitting a square field operator T, see 
Bouleau [2]. We denote VA the domain of the associated generator. 

Theorem 3. Let (Q,A, P, D, T) fre an error structure. Let be f G O fl L°° sncn t/iot 
/ > 0, E/ = 1. H^e pni P x = /.P. 

aj T/ie bilinear form £ 1 defined on VA n L°° fry 



/«A[«] + it;r[«,/] 



^[uju] = — E 
closable in L 2 (Pi) and satisfies for u,v G "DA fl L 



(5) 



5i[u,u] = -(Am,n) = -(n,An) = ^E[/r>,n]] (6) 

■wnere Ai[n] = A[n] + ^r[n, /]. 

6j Let (Di,£i) &e tae smallest closed extension of (VA fl L°°,£i). Then D C D 1; 
£i /oca/ and admits a square field operator T 1; and 

T 1 = T on D 

m addition VA C 2Mi and A x [u] = A[u] + ^r[n, /] /or a// n G VA. 

Proof. 1) First, using that the resolvent operators are bounded operators sending 
L°° into VA fl we see that VA fl L°° is dense in D (equipped with the usual 
norm (\\.\\ 2 L2 + £[]y/ 2 ), hence also dense in L 2 (F 1 ). 

2) Using that D fl L°° is an algebra, for n, v G PA fl L°° we have 

£x[n,n] = -E[fvA[u] + \vY[uJ}} = ^E[r[/n, u] - vF[u, /]] = ^E[/r>,n]]. 

So, defining A x as in the statement, we have Vn, n G "DA fl L°° 

5i[u,u] = -EifuAiu] = -EifuAiu]. 
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The operator A\ is therefore symmetric on VA n L°° under Pi. Hence the form E\ 
defined on VA fl L°° is closable, (Fukushima et al. [5], condition 1.1.3 p 4). 

3) Let (H>i,£i) be the smallest closed extension of (VA fl L°°,£i). Let be u G D 
and G PA fl with u n — > u in D. Using £i[/u n — « m ] < ||/||oo^[ M n — J* m ] an d 
the closedness of £\ we get u n — > u in D 1; hence D C D x . Now by usual inequali- 
ties we see that T[u n ] is a Cauchy sequence in L 1 (Pi) and that the limit Ti[u] does 
not depend on the particular sequence (u n ) satisfying the above condition. Then 
following Bouleau [2], Chap. Ill §2.5 p. 38, the functional calculus extends to Di, 
the axioms of error structures are fulfilled for (ft, A, Pi, Di, r^) and this gives with 
usual arguments the b) of the statement. □ 

1.3. Rajchman measures. 

In the whole paper, if x is a real number, [x] denotes the entire part of x and 
{x} — x — [x] the fractional part. 

Definition 1. A measure ji on the torus T 1 is said to be Rajchman if 



The set of Rajchman measures 7Z is a band : if \i G 7Z and if v <^ then v G 1Z, 
cf. Rajchman [18] [19], Lyons [15]. 

Lemma. Let X be a real random variable and let ^x(u) = Ke tuX be its character- 
istic function. Then 




when \n\ f oo. 





Proof, a) If lim| uKoo V x (u) = then ^ x (2vrn) = (F {x} }(n) ->■ 0. 
b) Let p be a probability measure on T 1 s.t. p G 71. From 




oo 



we have 




with a p (n) = 
x p p(dx) G 71, so 



f x p e 
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Now if P{x} e 7£, since l{ xe [^ p+1 [ } .F{ X } < we have 

lim E[e 2 ^]= lim £ E[e 2 ^l { x e[p , P+ i [} ] 

P 

which goes to zero by dominated convergence. □ 

A probability measure on R satisfying the conditions of the lemma will be called 
Rajchman. 

Examples. Thanks to the Riemann-Lebesgue lemma, absolutely continuous mea- 
sures are in 7Z. It follows from the lemma that if a measure v satisfies v*- ■ -*v e 7Z 
then u E 71. There are singular Rajchman measures, cf. Kahane and Salem [13]. 

The preceding definitions and properties extend to T d : a measure \i on T d is 
said to be in 1Z if fx{k) — > as k — > oo in 7h d . The set of measures in 1Z is a 
band. If X is Revalued, limiui^oo Ee^ u,x ^ = is equivalent to G 7£ where 

{x} = ({x^,. . .,{x d }). 

II. Finite dimensional cases. 

In the whole article =4> denotes the convergence in law, i.e. the convergence 
of the probability laws on bounded continuous functions. The arbitrary functions 
principle may be stated as follows: 



Proposition 1. Let X, Y, Z be random variables with values in R, R, and R m resp. 
Then 

({nX + Y},X,Y,Z) =4 (U,X,Y,Z) (7) 

where U is uniform on the unit interval independent of (X, Y, Z), if and only if fx 
is Rajchman. 

Proof. If /i is a probability measure on T 1 x R m , let us put 




2lnkx+{c ' y) fi{dx,dy), 



then fi n =^tiffi jl n (k, () ->■ p,(k, () VA; G Z, V( G R m . 
a) If Fx G 72. 

P({nX + y } ,x,y,z)(A;, Ci, C2, Cs) = E[exp{2*vrA;(nX + F) + + < 2 Y + i(C 3 , £)}] 

= J e 2i * knx f(x)F {x} (dx) 

with /(x) = E[exp{2i7rA;r + i&X + i( 2 Y + i(( 3 ,Z)}\{X} = x\. The fact that 
f-F{x} G 7?. gives the result. 



6 



b) Conversely, taking (£,(1,(2,(3) = (1,0, -2n, 0) gives P{x}(n) ->■ i.e. P x G 



Let us suppose now that Y is an K. -valued random variable, measured with 
an equidistant graduation corresponding to an orthonormal rectilinear coordinate 
system, and estimated to the nearest graduation component by component. Thus 
we put 

Y n = Y + -6{nY) 
n 

with 6{y) — {\ — {yi}, ■ ■ ■ , \ — {yd})- Let us emphasize that Y n is a deterministic 
function of Y. 

Theorem 4. a) Iffy is Rajchman and if X is W 71 -valued 

(X,n(Y n -Y)) =4 {X,(V u ...,V d )) (8) 

where the Vi 's are independent identically distributed uniformly distributed on (— |, |) 
and independent of X . 
For all up G C 1 n 

=4 (x^v^On) (9) 

i=l 

n 2 E[(v,(r n ) - v?0n) 2 |V = y] ^ ^ E Ll (^) ( 10 ) 

i=l 

m particular 

n*EMY n ) - V {Y)f] -+ E y [i^f(y)]. (11) 

i=i 

&,) zs of class C 2 , the conditional expectation n 2 W\<p{Y n ) — (p(Y)\Y = y] 
possesses a version n 2 (uo(y + -6{ny)) — <f(y)) independent of the probability measure 
P which converges in the sense of distributions to the function ^ A (p. 

c) //Py < dy on R d , G ^([0, 1]) 

(x,V(n(y n -y))) =4(x,v>G0). (12) 

d) We consider the bias operators on the algebra C 2 of bounded functions with 
bounded derivatives up to order 2 with the sequence a n = n 2 . IfPy&lZ and if one 
of the following condition is fulfilled 

i) \/i = l,...,d the partial derivative <9jPy in the sense of distributions is a 
measure <C Py of the form p^Py with pi G L 2 (Py) 7 

ii) Py = hl G fy with G open set, h G H 1 n L°°(G), h>0, 
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then hypotheses (HI) to (H4) are satisfied and 



= Ta A V + k E V?;P; case 
=^ A V^ + ^iE^yi case 



Proof. The argument for relation (8) is similar to the one dimensional case stated 
in proposition 1. The relation (9) comes from the Taylor expansion tp(Y n ) — tp(Y) = 

= EiLiC^M ~ Yi) Jo fi0^n,i, ■ ■ ■ , Y n ,i-i, Yi + t(Y n j — Yi), Y i+ i, . . . , Yd) dt 
and the convergence 



thanks to (8) and the following approximation in L l 



E 



V 9(nYj<fil(Y) - V 9{nY^) [\' t (. . . ,Y t + t(Y n , — Yi),. . .)dt 

Jo 



->■ 0. 



To prove the formulas (10) and (11) let us remark that 
n 2 E[^(Y n ) - V {Yf\Y = y] = 

= E V OinYi) f y^(. ..,Y l + t(Y n , -Y^,.. .)dt 
Jo 

y^Oinyi) / <Pi(yi + -0(nyi), ■ ■ + t-0(nyi), . . .)dt 
t Jo n n 



Py-a.s. 



each term (9(nyi) </^(. . .)dt) 2 converges to J 9 2 (pf(y) = ^y?f in L 1 and each term 
6(nyi)6(nyj) Jq ■ ■ ■ Jq ■ ■ ■ goes to zero in L 1 what proves the part a) of the statement. 

The part b) is obtained following the same lines with a Taylor expansion up to 
second order and an integration by part thanks to the fact that ip is now supposed 
to be C 2 . 

In order to prove c) let us suppose first that Py = lr 0j i]d.dy. Considering a 
sequence of functions ipk € C& tending to ip in L 1 we have the bound 

^ e i(u,X) e ivtl>(e(nY))j _ ^ e i(u,X) e roy> fc (0(nY))]| 

< \v\ J \^(9(ny)) - M9(ny))\dy 

= M E^ri • • • Ip' +1 ■ • • \1>(o( n yi) ■■■)- M°( n yi) ■ ■ -)\dyi ...dy d 
= HE---E/---/W(xi),...)-^(x 1 ),...)l^---^ 1 

= \v\\\ij) - tpk\W- 
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This yields (12) in this case. Now if P y < dy then P {y} < dy on [0, l] d and the 
weak convergence under dy on [0, l] d implies the weak convergence under f{y} what 
yields the result. 

In d) the point i) is proved by the approach already used in Bouleau [3] consisting 
of proving that hypothesis (H3) is fulfilled by displaying the operator A thanks to an 
integration by parts. The point ii) is an application of the Girsanov-type theorem 
3. □ 



Remarks. 1) About the relations (9) (10) (11), let us note that with respect to the 
form 

i 

when it is closable, the random variable £V appears to be a gradient : if we 
put ip* = £V then a we have 

i 

the square field operator associated to £. We will find this phenomenon again on 
the Wiener space. 

2) Approximation to the nearest graduation, by excess, or by default. When the 
approximation is done to the nearest graduation, on the algebra C 2 the four bias 
operators are obtained in theorem 4 with the sequence a n = n 2 , (with a n = n the 
four bias operators would be zero). 

We would obtain a quite different result with an approximation by default or by 
excess because of the dominating effect of the shift. 

If the random variable Y is approximated by default by Y^ = then 

n(Y^ -Y)=U-U and E[n{Y^ - Y)] -)• -- 

as soon as Y is say bounded. With this approximation, if we do not erase the shift 
down proportional to — and if we take a n = n we obtain first order bias operators 

without diffusion : A[(p] = —\ip' = — A[ip] and A = 0. The same happens of course 
with the approximation by excess. 

3) Extension to more general graduations. Let Y be an Revalued random vari- 
able approximated by Y n = Y + £ n (Y) with a sequence a n t oo on the algebra 
V = C{e^ u ' x \ u E R d }, the function £ n satisfying 

' a n E[\U 3 (Y)] -+ 

a n E[ip(Y)(u,£ n (Y)} 2 } ->■ E Y [ip.u*ju] Vv? eV,Vue R d 
^ < with 7y e L°°(P Y ) and ^ in distributions sense G L 2 (F Y ) 

k a n E[<p(Y){u,UY))]^0 W^eV. 
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Under these hypotheses we have 

Theorem 4bis. a) (HI) is satisfied and A[tp] = § Ylij lij Qxtx • 

b ) If for % — 1, . . . , d, the partial derivative <9jPy in the sense of distributions is a 
bounded measure of the form p«Py with pi G L 2 (Py) then assumptions (HI) to (H4) 
are fulfilled and Wip G V 



Proof. The argument is simple thanks to the choice of the algebra V and consists 
of elementary Taylor expansions to prove the existence of the bias operators. Then 



Historical comment. 

In his intuitive version, the idea underlying the arbitrary functions method is 
ancient. The historian J. von Plato [16] dates it back to a book of J. von Kries [12]. 
We find indeed in this philosophical treatise the idea that if a roulette had equal 
and infinitely small black and white cases, then there would be an equal probability 
to fall on a case or on the neighbour one, hence by addition an equal probability to 
fall either on black or on white. But no precise proof was given. The idea remains 
at the common sense level. 

A mathematical argument for the fairness of the roulette and for the equi- 
distribution of other mechanical systems (little planets on the Zodiac) was proposed 
by H. Poincare in his course on probability published in 1912 ([17], Chap. VIII §92 
and especially §93). In present language, Poincare shows the weak convergence of 
tX + Y"mod 2-7T when t f oo to the uniform law on (0, 27r) when the pair (X, Y) 
has a density. He uses the characteristic functions. His proof supposes the density 
be C 1 with bounded derivative in order to perform an integration by parts, but 
the proof would extend to the general absolutely case if we were using instead the 
Riemann-Lebesgue lemma. 

The question is then developed without major changes by several authors, E. 
Borel [1] (case of continuous density), M. Frechet [4] (case of Riemann-integrable 
density), B. Hostinski [9] [10] (bidimensional case) and is tackled anew by E. Hopf 
[6] , [7] and [8] with the more general point of view of asymptotic behaviour of 
dissipative dynamical systems. Hopf has shown that these phenomena are related 
to mixing and belong to the framework of ergodic theory. 



III. Infinite dimensional extensions of the arbitrary functions prin- 
ciple. 




the square field operator is F[ip\ = Ylijlij 



dip dip 
dxi dxj ' 



theorem 1 applies. 



□ 
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III.l. Rajchman type martingales. 

Let (Tt) be a right continuous filtration on A, P) and M be a continuous local 
{Tti P)-martingale nought at zero. M will be said to be Rajchman if the measure 
d(M, M) s restricted to compact intervals belongs to 1Z almost surely. We will show 
that the method followed by Rootzen [21] extends to Rajchman martingales and 
provides the following 

Theorem 5. Let M be a continuous local martingale which is Rajchman and s.t. 
(M, M)oo = oo. 

Let f be a bounded Riemann-integrable periodic function with unit period on H. 
s.t. J f(s)ds = 0. Then for any random variable X 

(X, [ f(ns)dM s ) ^ (X,W limM , M) ), (13) 
Jo 

the weak convergence is understood onM, x C([0, 1]) and W is an independent stan- 
dard Brownian motion. 

Before proving the theorem, let us remark that it shows that the random measure 
dM s behaves in some sense like a Rajchman measure. Indeed if Py G 1Z we have 

/y r 1 ry 

g(nx)F Y {dx) ^ g(x)dx F Y (dx) 
-oo Jo J —oo 

as soon as g is periodic with unit period, Riemann-integrable and bounded. Now 
applying the theorem to the Brownian motion gives the similar relation 

f f(n8) dB s =4 ( f f 2 (s)d S y/ 2 f dW s . 
Jo Jo Jo 

Proof. We consider the local martingale N t = J ' f{ns)dM s . 

a) In order to be sure that (N, N)^ = oo, we change N t into N t 
with f n (s) = f(ns) for s e [0,1), f n (s) = for s e [l,n] and f n (s) 
We put S n (t) = inf{s : (N,N) S > t}. 

b) We want to show 

nmN Sn )\ -+ nmw)\ v^gl^p) vfgc^i]). (u) 

It is enough to consider the case £ > 0, E£ = 1, and £ may be supposed to be 
J-^-measurable for a deterministic time T large enough. Let be P = £.P and D(t) = 
E^IJ^]. The process 

Mt = M t - [ D-\s)d{M,D c ) s 
Jo 



= f*f n (s)dM s 
— 1 for t > n. 
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~ S (t) ~ 

is a continuous local martingale under P. Therefore J Q f n (s) dM s is a Brownian 
motion under P (Revuz and Yor [20] p. 313 theorem 1.4 and p 173). Writing 



/ f n {s)dM s = f n (s)dM s + ^-> d (M,D c ) s 

Jo Jo Jo ^K 8 ) 

and noting that d{M,D c ) s vanishes on ]T, oof, in order to show (14) it suffices to 
show 

** US) 



-S n (t) 



hence to show 



sup 

0<t<T 



sup 

0<t<l 



D(s) 

1 f(ns) 



d(M, D c 



d{M, D c 



a.s. when n — > oo 



— > a.s. when n — > oo 



/o D(s) 

and, because M is Rajchman this comes from the following lemma : 
Lemma. Let f be as in the statement of the theorem, then V/x G 1Z 



sup 

0<«1 



/ f(ns)fi(ds) 
Jo 



— > as n — > oo. 



Proof. We have 



f(ns)fx(ds) ->■ / f(s)ds / //(ds) = 0. 
Jo Jo 



Since / is bounded, the functions J* f(ns)fj,(ds) are equi-continuous and the result 
follows from Ascoli theorem. □ 

c) This proves the following stable convergence 

(X, / f(ns)dM s ) ^ {X,W) 
Jo 

and by the fact that the limit 

,t ,i 
/ 2 (ns)d(M,M) s ^ / f 2 (s)ds(M,M} t 
Jo 

is a continuous process, this gives the announced result. □ 
Remark. If f(s)ds ^ 0, then keeping the other hypotheses unchanged, we obtain 



V. / f(ns)dM s ) =4 (x,( ^ f(s)ds)M. + ( f\f - f f) 2 ) 1/2 W {M , 
'o V Jo Jo Jo 



M). 
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We study now the induced limit quadratic form when the martingale M is ap- 
proximated by the martingale M™ = M t + J* ^f(ns)dM s . The notation is the same 
as in the preceding section and / satisfies the same hypotheses as in theorem 5. 

Theorem 6. Let M be a Rajchman martingale s.t. M\ e L 2 and rj, ( bounded 
adapted processes. Then 



n 2 E 



(exp{* I rjsdM?} - exp{i I Vs dM s })(exp{i [ ( s dM?} - exp{i I ( s dM s }) 
Jo Jo Jo Jo 



->■ — E 



exp{i / (r) s + (s)dM s } [ r) s ( s d(M,M) t 
Jo Jo 



f(s)ds. 







Proof. By the fundamental formula of calculus (finite increments formula), the first 
term in the statement may be written 

-E[exp{* l (r) s + QdM a } [ r] s f(ns)dM s [ ( s f(ns)dM s ] + o(l) 
Jo Jo Jo 

therefore, thanks to theorem 5, the statement is a consequence of the following 
lemma : 

Lemma. Suppose EM 2 < oo and i] adapted and bounded, then the random variables 
r) s f(ns)dM s are uniformly integrable. 

Proof. It suffices to remark that their L -norm is equal to E f* r] 2 J 2 {ns) d(M, M) s 
hence uniformly bounded. □ 

III. 2. Sufficient closability conditions on the Wiener space. 

The closability problem of the limit quadratic forms obtained in the preceding 
section, may be tackled with the tools available on the Wiener space. 

Let us approximate the Brownian motion (S t ) t£ [ 0il ] by the process = B t + 
Jo nf( ns ) dB s where / satisfies the same hypotheses as before. 

Theorem 7. a) Let £ e L 2 ([0, 1]), and let X be a random variable defined on the 
Wiener space, i.e. a Wiener functional, then 



X,n(e W {ij\dB n }-exp{ij\dB})^j =4 (x,\\f\\ L2 (exp{i j\dB}) 



# 



(15) 

here for any regular Wiener functional Z we put Z # (u, w) = J D s ZdW s , where W 
is an independent Brownian motion. 
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b) 

n 2 E [ (e iC.B» _ e ie.B )2 ] ^ _ E [ e ^.B] f 1 ^ds\\f\\ 2 L2 (16) 

on £/ie algebra £{e^ B } the quadratic form — ^'E[e 2l ^' B ] £, 2 ds is closable, its closure 
is the Ornstein- Uhlenbeck form. 

Proof, a) The first assertion comes easily from the similar result concerning Ra- 
jchman martingales using the fact that e ia h n^ ns ^ dBs da — > 1 in L' p p e [1, oo[. 

b) The obtained quadratic form is immediately recognized as the Ornstein- 
Uhlenbeck form which is closed. It follows that hypothesis (H3) is fulfilled and 
the symmetric bias operator is 

A[e^n = (-*- JtdB-\j e'SW □ 

If instead of the Wiener measure m, we consider the measure mi = h.m for 
an h > 0, h G D ou H L°° where D ou (= D 2,1 ) denotes the domain of the Ornstein- 
Uhlenbeck form, we know by the Girsanov-type theorem 3 that the form 



is closable, admits the same square field operator on D ou , and that its generator A\ 
satisfies 

A 1 [<p] = A[<p] + — V ou [if, h] for <p e VA ou 

Since the point a) of the theorem is still valid under mi because of the properties of 
stable convergence, the preceding theorem is still valid under mi, the Dirichlet form 
being now 

£i[ifi] = -EiF^]] for V eB ou . 

Remark. Let us come back to the general case of Rajchman martingales. If we 
suppose the Rajchman local martingale M is in addition Gaussian, which is equiva- 
lent to suppose (M, M) deterministic, then on the algebra C{e l ^ dM ; £ deterministic 
bounded } the limit quadratic form 



_ E[e iJ(v+OdM (\ sr}sd{M ,M) s \\\f\\h 
Jo 



is closable, hence (H3) is satisfied. 

Indeed, it suffices to exhibit the corresponding symmetric bias operator. But 
by the use of the calculus for Gaussian variables, it is easily seen that the operator 
defined by 

Apt*™] = jS*™ JtdM-\J ed(M,M))j J f 2 ds 
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satisfies the required condition. 



□ 



III. 3. Approximation of the Ornstein-Uhlenbeck gradient. 

Let m be the Wiener measure on C([0, 1], R). Let be a real periodic function of 
period 1 such that 9(s)ds = et 9 2 (s)ds = 1. We consider the transformation 
R n of the space L^(m) defined by its action on the chaos: 
if X = | si< ... <sfc f( Sl , . . . , s k )dB Sl . . . dB Sk for / G L% m ([0, l] k , C), 

R n (X) = [ f( Sl , s k y^ ns ^dB Sl . . . e^ ns ^dB Sk . 

Since \\X\\ 2 L2(m) = j si< ^ <Sk \f\ 2 d Sl ...ds k = £11/11^, R n is an isometry from 
L&m) into itself and V^4([0,l]) 

|| e /^-|/^ s || i2 = ehfW 2 ^. 



From the relation 



k „i 

n{e iT.U d ^) - 1) = iJ20(ns p ) / e<^r> e{ns ^da 

p=i Jo 



it follows that if X belongs to k-th chaos 

\\n(R n (X) - X)\\l 2 < k^XfMl 
then, denoting A the Ornstein-Uhlenbeck operator, for X e f(^4) 

||n(i?„(X)-X)|| L2 < 211^X1111^1100 

and we can state: 

Theorem 8. If X e V(A) 

(-in(R n (X)-X),B) =4 (X#,5) 

mtt X # = / X D s XdW s where W is an independent Brownian motion. 

Proof. If X belongs to the /c-th chaos, expanding the exponential by its Taylor 
series gives 

n(R n (X) -X)=i f ...,s k )J2 0(ns p )dB Sl . . . dB Sk + Q n 

Js 1 <-<s k p=1 
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with \\Q n \\ 2 <±k 2 \\e\a\X\\ 2 . 

Now, since j Sl< ... <Sp< ... <Sk H s i> • • • > s k )6{ns p )dB Sl . . . dB Sp . . . dB Sk converges sta- 
bly to J Sl< ... <Sp< ... <Sk h(si, . . . , s k )dB Sl . . . dW Sp . . . dB Sk we obtain that 

It< S2 <-< Sk /(*> s 2, • • • , s k )dW t dB S2 . . . dB Sk 
-in(R n (X) L > > - /(Sl ' S ^ dB - dWt ■ ■ ■ dB °* 

+ f sl< ... <Sk _ 1<t /(*!> • • • > s fc-i. O^i • • • dB Sk _ x dW t 

which is equal to / D s (X)dW s = X* . 

For the general case, we approximate X by X k for the norm D 2,2 and reasoning 
with the characteristic functions yields the result (see the proof of theorem 10 be- 
low). □ 

By the properties of the stable convergence, the convergence in law of theorem 
8 still holds under fh <C m. 

Theorem 9. VX e V(A) 

n 2 E[\R n (X) - X\ 2 } -> 2S[X] 
where £ is the Dirichlet form associated with the Ornstein- Uhlenbeck operator. 

Proof. As R n preserves the chaos and the expansion of n(R n (X) —X) on the chaos 
is dominated by that of 2||0|| OO AX', it suffices to argue when X is in the fc-th chaos. 
Starting from 



n 2 E\R 



t (X) -X\ 2 = n 2 [ |/| 2 | e * -l\ 2 d Sl ... ds k 



' Sl<---<S k 

expanding the exponential and estimating the remainder we obtain 

\imn 2 E\R n (X) - X\ 2 = [ \f\ 2 d Sl ...ds k [ 9 2 dt = k\\X\\ 2 

k ' J[o,i] k Jo 

what gives the result. □ 

Following the same arguments, it is possible to show that the theoretical and 
practical bias operators A and A defined on the algebra C{e^ dB ; £ G C 1 } by 

n 2 E[(R n (X) -X)Y] = (AX,Y) L 2 {m) 
n 2 E[(X - R n {X))R n {Y)\ = (AX,Y) LHm) 

exist and are equal to A. 
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III. 4. Isometries on the Wiener space. 

Let us now consider a d-dimensional Brownian motion (B t ). 

Let 1 1 — y M t be a deterministic bounded measurable periodic map with period 1 
with values in the space of d x d othogonal matrices such that j Q M s ds = 0, (for 
instance a rotation of angle 2nt). We denote still m the Wiener measure law of 
(B t ). The transformation B t h-> J* * M s dB s induces an endomorphism T M isometric 
in L p (m), 1 < p < oo. We put M n (s) = M ns and T n = T Mn - 

Theorem 10. Let X be in L l {m). Let be m Cm, we have under m: 

(T n (X),B) =4 

The convergence in law is understood on R x C([0, 1]) and X(u>) denotes a random 
variable with the same law as that of X under m, function of a Brownian motion 
W independent of B. 



Proof, a) If X has the form 

X = exp{* [\.dB + l ffcfds} 
Jo 1 Jo 

for an element £ G L 2 ([0, l],R d ), we have 

T n (X)=exp{iJ £M n (s)dB a + ±J \£\ 2 ds}. 

where £* denotes the transposed of £ s . If we put Z r t l = j^ t s M n {s)dB s 

(Z n ,Z n ) t = f e s M n (s)M* n ( S )t s ds = fltf^ds 
Jo Jo 

is a continuous function. Now by theorem 4 

C s M n (s)ds^ [ C s ds [ M n (s)ds = 0. 



Since the functions t \-t f* £*M n (s)ds are uniformly continuous (M is bounded) 
by Ascoli theorem sup t \ f* £*M n (s)ds\ — > 0. The argument of Rootzen applies once 



more 



(f'CM n dB,B) =4 {f'i.dW^B) 
Jo Jo 



giving the result by the continuity of the exponential function. 

b) For X e L 1 {m), we consider X k linear combination of exponentials of the 
above form approximating X in L l {m). 
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By a) we have \/h e L 2 {[0, l],R d ) 

but 

m^T n (X) e ifh.dB _ E [ e iuT n (X k ) e ifh.dB^ < \ u \E\T n (X)-T n (X k )\ = \u\ \\X-X k \\ L i 

what gives the result. 

c) This extends to m < m by the properties of stable convergence □ 

III. 5. Stochastic differential equations from dynamics. 

In the case f(x) = 9{x) — \ — {x}, the approximation used in parts III. 1 and 
III. 2 approaches B t by 

ft 



and yields limit of the type 



n(B. - B% B.) = („ /""(a - J- - ^)dB s , B) (^=W„ B) 

Jo 2n n vl2 



Hud n\ ( 1 

ynyo. - id ): idj = yn - — - 

and 

(n j\ 8 - [ ^)dB s , n J\B a - B^)ds, B?j =4 (JL^.+Ifl , - JLwr+Ifl , B) 

(18) 

Now, when we solve by the Euler method a stochastic differential equation of the 
type defining a diffusion process and expand the coefficients in series, we encounter 
integrals of the type 

[\ s -M) d B a , ['(B s -B^)ds 
Jo n Jo 

but also of the type 



(B s -Bw)dB s (19) 

n 

and these last ones, by the central limit theorem, yield the convergence 

(v^J o (B s -B^)dB s ,B}j =U(±W.,,B.). (20) 

Then let us remark that 

a) the limits (18) are generally hidden by the limits (20) because of the order of 
magnitude of the coefficients n and \Jn respectively, 

b) in (20), the conditional law of the random variable f Q B\ ns } dB s with values 

n 

in C([0, 1],M) given f Q B s dB s is not reduced to a Dirac mass, the approximation is 
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not deterministic (as can be seen, for instance, by changing the sign of a Brownian 
path after the time T = infls) 2 ^-, B s = 0}, the a-field generated by f B s dB s being 
°~{B 2 ,s < 1) ). In (17) instead, B n is a deterministic function of B. 

Nevertheless, for some stochastic differential equations the limits (18) remain 
dominant and determine the convergence. This concerns, for instance, stochastic 
differential equations of the form 

X\ = xl + ft f n (X 2 )dB s + ft / 12 (Xj, Xl)ds 

where I is w ith values in R k \ X 2 in M fe2 , B in R d and f ij are matrices with 
suitable dimensions. Such equations are encountered to describe the movement of 
mechanical systems under the action of forces with a random noise, when the noisy 
forces depend on the position of the system and the time. Typically 

X t = X + fiv a ds 

V t = V + f a(X s , V s , s)ds + f* b(X s , s)dB s 

which is a perturbation of the equation = a(x,^,t). In such equations the 
stochastic integral may be understood as Ito as well as Stratonovitch. For the equa- 
tion (21) the iterative method of Kurtz and Protter [14] (see also Jacod and Protter 
[11]) may be applied starting with the results obtained in generalizing of the arbi- 
trary functions principle. This yields the following result that we state in the case 
ki = &2 = d = 1 for simplicity. 

Theorem 11. If functions f lj are C\, and if X n is the solution of (22) by the Euler 
scheme, 

(n(X n -X),X,B)^(U,X,B) 
where the process U is solution of the stochastic differential equation 

U ^ = E f d -^(X s )U k s dYi - W* §^(X S ) £ f m (X s )dZr 
kj Jo ° Xk kj J ° OXk 

where Y s = (B s , s)' and 

dZf = -^dW s + \dB K 



dZf = ~^dW s + \dB s 
dZf = f 

and as ever W is an independent Brownian motion. 

Thus the Euler scheme for solving this kind of equations encountered in me- 
chanics gives rise to an asymptotic weak limit, but in - and based on the arbitrary 
functions principle, instead of being in ^ and based on a version of the central limit 
theorem. 
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